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ABSTRACT 

We report the discovery of four kpc- scale binary AGNs. These objects were originally selected from the 
Sloan Digital Sky Survey based on double-peaked [Om] AA4959,5007 emission lines in their fiber spectra. 
The double peaks could result from pairing active supermassive black holes (SMBHs) in a galaxy merger, 
or could be due to bulk motions of narrow-line region gas around a single SMBH. Deep near-infrared (NIR) 
images and optical slit spectra obtained from the Magellan 6.5 m and the APO 3.5 m telescopes strongly 
support the binary SMBH scenario for the four objects. In each system, the NIR images reveal tidal features 
and double stellar components with a projected separation of several kpc, while optical slit spectra show two 
Seyfert 2 nuclei spatially coincident with the stellar components, with line-of-sight velocity offsets of a few 
hundred km s" 1 . These objects were drawn from a sample of only 43 objects, demonstrating the efficiency of 
this technique to find kpc- scale binary AGNs. 

Subject headings: black hole physics - galaxies: active - galaxies: interactions - galaxies: nuclei - galaxies: 
Seyfert - quasars: general 



1. INTRODUCTION 

Despite decades of searching, and strong theoretical rea- 
sons to believe they exist, merging active supermassive black 
holes (SMBHs) remain difficult to find. Here we present the 
first four secure candidates from our ongoing search for kpc- 
scale binary 5 AGNs. 

It is widely apprec iated that galaxies are built up hierar- 
chically via mergers (iToomre & Toomrelll972|) . Since most 
massive galaxies are believed to harbor a central SMBH 
(Kormendy & RichstoneJ Il995|) . galaxy merg ers would re- 
sult in the formation of bin ary SM BHs dBegelman et al.l 
1980; Milosavljevic & Merritt 2001). Major mergers be- 
tween galaxies are responsible, at least in some mod- 
els, for triggering black hole accretion activity in the 
most luminous AGNs. Merger-based models do a de- 
cent job of fitting observational data on local SMBH de- 
mogr aphics and AGN statisti cs at different cosmic epochs 
(e.g., iKauffmann & HaehneE 120001 : IVolonteri et aD 120031 : 
IWvithe&Loebll2003l: iHopkins et al.1 120081: [Shen 2009) , and 
on th e core properties of elliptical galaxies (e.g. jFaber et al.l 
11997b iKormendv & Bended l2009h . but direct observational 
evidence for binary SMBHs remains s urpris i ngly scarce (e.g., 
iKomossa et all 120031: fBianchi et al. 200 |; iRodriguez et al.1 
120061 ICom erford et al. 2009b; G reen et alJl2010h. 

iGerke et al.1 (l2007h and IComerford et al. ( 2009a) reported 
two binary-A GN candid ates in the DEEP2 galaxy redshift 
survey dDavis et al .1120031) . which show spatially resolved (on 
~ kpc scales) double-peaked [Om] A5007 emission with 
velocity offsets of a few hundred km s -1 . We have selected 

1 This paper includes data gathered with the 6.5 meter Magellan Tele- 
scopes located at Las Campanas Observatory, Chile, and with the Apache 
Point Observatory 3.5-meter telescope, which is owned and operated by the 
Astrophysical Research Consortium. 

2 Department of Astrophysical Sciences, Princeton University, Peyton Hall 
- Ivy Lane, Princeton, NJ 08544 

3 Princeton-Carnegie Fellow 

4 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS-51, 
Cambridge, MA 02138 

5 We refer to both a bound pair of SMBHs in a Keplerian orbit, and a 
pair of SMBHs at large separations where the galactic potential dominates, 
as binaries. 



a similar sample of 167 do uble-p eak ed TO ml na rrow-line 
AGNs (see also ISmith et ai1l2009l and IWang et al.ll2009l for 
similar statistical studies, as well as work on individual 
interesting objects, e.g. Xu & Komossal l2009h from the 
Seventh Data Release of the Sloan Digital Sky Survey 
(SDSS; lYorketalJl200a lAbazajian et aDl2009f) and studied 
their host-gala xy and em ission-line properties and selection 
completeness dLiu et al.ll2010|) . The double-peaked features 



such as outflows or rotating disks (e.g., 


Axon et al.l 11998b 


Veilleux et al.1 1200 ll: ICrenshaw et al.1 12010 


), or to a pair of 



AGNs. To discriminate between these scenarios for individual 
objects and to identify bona fide AGN pairs, we are conduct- 
ing deep near-infrared (NIR) imaging and spatially resolved 
optical spectroscopy for the objects in the double-peaked 
narrow-line sample. Here we report initial results on the dis- 
covery of four binary AGNs, SDSSJ1 10851.04 + 065901.4, 
SDSSJ1 13 126.08 - 020459.2, SDSSJ1 14642.47 + 5 1 1029.6, 
and SDSSJ133226.34 + 060627.4. Deep NIR imaging re- 
vealed double stellar bulges in SDSSJ1 1085 1.04 + 065901. 4, 
SDSSJ1 13126.08 - 020459.2, and SDSSJ133226.34 + 
060627.4, while SDSSJ1 14642.47 + 51 1029.6's SDSS image 
shows clear evidence of a recent merger. We have obtained 
slit spectroscopy of all four sources . We describe our obser- 
vations data analysis in § 32.1112.21 follo wed b y discussion 
of the nature of the ionizing sources in 32.31 We discuss 
our results and conclude in ^3j A cosmology with £l m = 0.3, 
ft A = 0.7, and h = 0.7 is assumed throughout. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. NIR Imaging 

We obtained K s and / images for SDSSJ1 13126.08 - 
020459.2, a K s image for SDSSJ1 10851.04 + 065901.4, 
and a / image for SDSSJ133226.34 + 060627.4 on the 
nights of 2009 December 29 through 2010 January 2 UT 
using the Per s son's A uxiliary Nasmyth Infrared Camera 
(PANIC; iMartini et al.1 120041) on the 6.5 m Magellan I 
(Baade) telescope. PANIC has a 2' x 2' field of view 
(FOV) and 0".126 pixels. The observing conditions were 
clear but not photometric, with seeing ranging between 
0".5 and 0".8 in the optical. Exposure times were 1080 
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Fig. 1.— (Left) NIR images for SDSSJ1 14642.47 + 511029.6 (K s ) 
from 2MASS and for SDSSJ110851.04 + 065901.4 (K s ), SDSSJ113126.08- 
020459.2 (Ks), and SDSSJ133226.34 + 060627.4 (J) from our PANIC obser- 
vations. North is up and east is to the left. Spatial scales are indicated on 
each panel. Colors are an asinh scale to intensity with a stretch to highlight 
the central nuclei. (Right) SDSS r-band images to show tidal features and 
galaxy-scale environment. 

s and 2160 s in K s for SDSSJ1 10851.04 + 065901.4 and 
SDSSJ1 13126.08 -020459.2 respectively, and 1080 s in / 
for both SDSSJ113126.08-020459.2and SDSSJ133226.34+ 
060627.4. We reduce d PANIC data with the Carnegie Super- 
nova Project pipeline (Hamuy et al. 2006) following standard 
procedures. 

Figure Q] shows the K s image for SDS SJ1 14642.47 + 
5 1 1029.6 from 2MASS (ISkrutskie et al.ll2006l) . and the K s im- 
ages for SDSSJ1 1085 1 .04+065901 .4 and SDSSJ1 13 126.08 - 
020459.2 and the / image for SDSSJ133226.34 + 060627.4. 
Also presented are the SDSS r-band images to illustrate the 
galaxy- scale environment and tidal features. The NIR images 
show two stellar bulges in each of these systems. In addition, 
all four systems show tidal features suggestive of mergers. For 
SDSSJ1 14642.47 + 51 1029.6, this is clearly seen in its SDSS 



image; in SDSSJ110851.04 + 065901.4, SDSSJ113126.08- 
020459.2, and SDSSJ133226.34 + 060627.4, the SDSS im- 
ages show tentative tidal features which are confirmed in deep 
PANIC images. The PANIC image of SDSSJ1 10851.04 + 
065901 .4 also reveals a small companion (~3 // west to the nu- 
clei) which is barely resolved in its SDSS image. The two stel- 
lar nuclei in SDSSJ1 1085 1.04 + 065901. 4 appear marginally 
resolved in the center; their SDSS counterparts exhibit differ- 
ent colors which could indicate different stellar populations or 
dust effects. We performed aperture photometry for the two 
stellar nuclei to infer their luminosity ratio. Table Q] lists the 
estimated luminosity ratios in K s or in /. Assuming similar 
mass-to-light ratios in the two stellar nuclei, they all appear 
to be major mergers with comparable stellar masses. Also 
listed in Table [T] are the projected separation A0;, the stellar 
mass of the whole galaxy es timated based on SDSS photom- 
etry dKauffmann et al.ll2003h . and the 2MASS K s magnitudes 
and ^-/colors. 

2.2. Optical Slit Spectroscopy 

Because of the lack of spatial information in the SDSS 
spectra, it is unclear whether the two velocity components 
seen in [Om] are spatially associated with the two stellar 
bulges. To resolve this question we obtained slit spectra for 
SDSSJ1 1085 1 .04 + 065901 .4, SDSSJ1 13 126.08 - 020459.2, 
and SDSSJ133226.34 + 060627.4 on the nights of 2010 Jan- 
uary 12 through 14 UT using the Low-Dispersion Survey 
Spectrograph (LDSS3) on the 6.5 m Magellan II (Clay) 
telescope. The observing conditions were clear but not 
photometric, with seeing ranging between 0".7 and I'M. 
LDSS3 has a 8'.3 diameter FOV and 0".188 pixels. We 
employed a l // x4 / long-slit with the VPH-Blue grism to 
cover H/3 and [Om] AA4959,5007, and the VPH-Red grism 
(with the OG590 filter) to cover Ha and [Nil]. We ori- 
ented the slit to go through the two stellar bulges seen 
in the NIR images (with PA = 141°, 80°, and 17° for 
SDSSJ1 1085 1 .04 + 065901 .4, SDSSJ1 13 126.08 - 020459.2, 
and SDSSJ133226.34 + 060627.4 respectively). The spec- 
tral resolution was 3.1 (6.3) A FWHM in the blue (red). 
Total exposure times were 4800s, 6300s, and 4500s for 
SDSSJ1 1085 1 .04 + 065901 .4, SDSSJ1 13 126.08 - 020459.2, 
and SDSSJ133226.34+060627.4respectively in the blue, and 
1800s for all three objects in the red. We took wavelength 
calibration spectra and flat fields after observing each object. 
We obtained spectra for SDSSJ1 14642.47 + 51 1029.6 on the 
night of 2010 February 15 UT using the Dual Imaging Spec- 
trograph (DIS) on the Apache Point Observatory 3.5 m tele- 
scope. The sky was photometric, albeit with poor seeing rang- 
ing between 1".7 and 2".7. DIS has a 4'x6' FOV and 0".414 
pixels. We adopted a 1".5 x 6' slit and the B 1200+R1200 grat- 
ings centered at 5500 and 7450 A. The slit was oriented with 
PA = 59° to go through both bulges shown in the SDSS image. 
The spectral resolution was 1.8 (1.3) A FWHM in the blue 
(red). The total exposure time was 4800s. We reduced the 
DIS and L PS S3 data following standard IRAF 6 procedures 
dTodvll 19861) and with the COSMOS reduction pipeline 7 . We 
observed white dwarfs EG21 and G191B2B for spectrophoto- 
metric calibration. We applied a telluric correction from our 

6 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 

7 http://www.ociw.edu/Code/cosmos 
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Fig. 2.— (Top four rows) LDSS3 and DIS 2D spectra over the H/3-[0 III] 
(left column) and Ho;-[N II] (right column) regions. Colors indicate flux den- 
sities with asinh scales and a stretch to illustrate both nuclei and extended 
emission. In both SDSSJ1 14642.47 + 511029.6 and SDSSJ133226.34 + 
060627.4, the blueshifted (redshifted) component corresponds to the SW 
(NE) stellar bulge in its image (Figure[[}. In SDSSJ1 13126.08 -020459.2, 
the blueshifted (redshifted) component corresponds to the E (W) bulge in 
its image, while in SDSSJ110851.04 + 065901.4, the blueshifted (redshifted) 
component corresponds to the NW (SE) bulge in its image. (Bottom row) 
Stellar-continuum subtracted spectra for SDSSJ133226.34 + 060627.4 from 
LDSS3 (red, dashed curves) and from SDSS (data in black, best-fit model in 
green, and the two velocity components in cyan and magenta). The vertical 
dotted lines are drawn at the systemic redshift of each object as measured 
from stellar absorption features. In SDSSJ133226.34 + 060627.4 the SDSS 
fiber covered only part of the redshifted component. 

standard stars after extracting ID spectra. 

The spatial correspondence between the two stellar con- 
tinuum peaks and two [Om] velocity peaks in all cases 
strongly indicates we are seeing two active galaxies (Fig- 
ure [2]). We see unambiguous spatial separations between the 
red and blue velocity components in all targets and in each 
strong emission line. We measure this spatial offset, A0 e , 
from a spatial profile summed over twice the FWHM of each 
line (although we combine Ha and [Nil] due to blending). 
We also checked these results by directly examining the 2D 
flux-density profile of each velocity component. The spa- 
tial offsets derived from each strong line [H/3, [O in] A4959, 
[Om] A5007 (Aftoni] in Table and the Ha-[Nn] region 
(A(9 H a in Table Q}] agree within ~0".3. Note in the case of 
SDSSJ1 14642.47 + 51 1029.6, the [O in] AA4959,5007 lines 
fall in the dichroic edge and are poorly detected in the red- 
shifted component 8 , but the spatial separation is clearly mea- 

8 In SDSSJ1 14642.47 + 51 1029.6, the blueshifted component itself shows 
double peaks which also appear spatially offset from one another; our mea- 
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Fig. 3. — Diagnostic line ratios measured from our LDSS3 and DIS spec- 
tra. In each system the two velocity components are plotted using the same 
color, with the filled (open) symbol denoting the blueshifted (redshifted) 
component. Gray scales indicate number densities of 31 ,179 emission- 
line galaxies from the SDSS DR4 (Kauffmann et al. 2003). Also shown 
is the theoretical starburst limit (dotted curve) from Kewley et al. ( 2001) 
along with the empirical separ ation (dashed curve) between H n regions and 
AGNs ( Kauffmann et al. 2003). Pure star-forming (SF) galaxies lie below 
the dashed curve, AGN-dominated objects lie above the dotted curve, and H 
n-AGN composites (C) lie in between. 

sured in the Ha and [Nil] lines. The separation between the 
two velocity components of the emission lines is consistent 
with that of the double stellar bulges (A0; in Table [T]). More 
directly, the emission lines are coincident with the stellar con- 
tinua seen in the spectra. This spatial coincidence suggests 
that the double-peaked emission line profiles in the SDSS 
spectra of these objects are due to the orbital motion of the 
gas associated with the double stellar bulges. This is further 
supported by the tidal features indicative of mergers seen in 
all four systems ( 32.11 Figure [T]). 

2.3. The Nature of the Ionizing Sources 

The nature of the ionizing source in each object is 
constrained using the BPT diagram dBaldwin et all fl98lt 
Veill eux & Osterbrocklll987|) featuring the emission line ra- 
tios [Om] A5007/H/3 and [Nil] A6584/Ha (Figure El). We 
fit multi-component models to the emission lines over the 
st ellar-continuum subtracted LDSS3 and DIS spectra (refer 
to lLiu et aTll2010L for details on constructing continuum mod- 
els and emission line fitting) and measured line ratios. For 
SDSSJ1 1085 1 .04 + 065901 .4, SDSSJ1 13 126.08 - 020459.2, 
and SDSSJ133226.34 + 060627.4, we extracted ID spectra 
summing over both velocity components as they are partially 
blended spatially, whereas for SDSSJ1 14642.47 + 51 1029.6 
we extracted ID spectra separately for each of the two ve- 
locity components. We compared measurements based on 
the SDSS and the LDSS3 and DIS spectra and found general 
agreement in the line ratios; the difference in line fluxes is 
likely due to the difference in the fiber and slit coverage. Fig- 
ure [2 (bottom panels) shows an example of our emission line 

surements for the blueshifted component represent the luminosity-weighted 
average of both peaks. The redshifted component overlaps with the redshifted 
peak of the blueshifted component in velocity space, and we extracted ID 
spectra for the two components separately to disentangle them properly. 
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TABLE 1 

Host-Galaxy and Emission-Line Properties. 



A0 f ASi A0 [Om] A6>Ha 
SDSS Designation Redshift M* K s g-r J-K s (") (kpc) L\/Li L [Qmh L [Qmh (") (") n ei n £l 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) 

SDSSJ1 10851.04 + 065901.4 0.1816 11.0 14.29 0.8 1.3 0.5 1.5 1.0 8.58 7.94 0.9 0.6 2A^\ 2.4+^ 

SDSSJ113126.08-020459.2 0.1463 11.3 14.51 1.0 1.3 0.6 1.5 1.8 7.81 7.72 0.6 0.6 2.7^ -0.3+^ 

SDSSJ1 14642.47 + 51 1029.6 0.1300 10.9 14.32 1.0 1.4 2.7 6.3 0.9 8.35 7.79 2.5 2.8 3.l4| 2.9^ 

SDSSJ133226.34 + 060627.4 0.2070 11.1 14.72 1.2 1.6 1.5 5.1 1.6 7.48 8.25 1.6 1.6 2.8^| 2.4™ : f 

Note. — The subscri pts "1" and "2" denote blueshifted and redshifted components. Col. 3: total stellar masses in the form of log(M* /M©) for the whole galaxy based on SDSS 
photometry from Kauffmann etal. ( 2003); Col. 4: from 2MASS; Col. 5: color from SDSS photometry; Col. 6: color from 2MASS; Cols. 7 & 8: projected separation based on NIR 
imaging; Col. 9: J (for SDSSJ133226. 34 + 060627. 4) or Ks-band (for the others) luminosity ratio of the two stellar bulges estimated with aperture photometry; Cols. 10 & 11: blueshifted 
and redshifted [O III] A5007 component luminosities in the form of log(L/L0 ) ; Col. 12: projected spatial offset measured from [O III] A5007; Col. 13: projected spatial offset measured 
from Hen and [Nil]; Cols. 14 & 15: electron density in the form of log(« e /cm~ 3 ) of the blueshifted and redshifted components. 



fits and of the difference in SDSS and slit spectra. In all cases, 
the diagnostic line ratios indicate Seyfert-type ionization 
(except for the redshifted component in SDSSJ1 14642.47 + 
51 1029.6, which could also be a LINER or a composite). Ta- 
ble[T]lists [Om] A5007 luminosities for both velocity compo- 
nents based on our LDSS3/DIS spectra. 

While the diagnostic line ratios indicate the presence of 
hard ionizing photons, it is ambiguous whether there are two 
ionizing sources. If only one nucleus were to ionize both gas 
components, the ionization parameter in the off-nuclear re- 
gion would be much lower than that in the nuclear region with 
similar electron densities, because the number density of ion- 
izing photons decays as r~ 2 . We measured electron densities 
n e using the diagnostic line ratio [S n] A6717/[S n] A6731. For 
SDSSJ110851.04 + 065901.4, SDSSJ1 14642.47 + 511029.6, 
and SDSSJ133226.34 + 060627.4, we found comparable n e 
values in the two velocity components (Table [T]). The compa- 
rable electron densities and high ionization states (Figure [3]) 
in both narrow line components suggest that there are sepa- 
rate hard ionizing sources associated with each line compo- 
nent. For SDSSJ1 13126.08 -020459.2, the inferred n e of 
the redshifted component is ~1000 times smaller than the 
blueshifted component, but the uncertainties of the measure- 
ment are so large that they may differ by no more than a fac- 
tor of ~ 10. Considering typical n e values of narrow-line re- 
gion gas in stellar bulges, in SDSSJ1 13126.08 -020459.2 the 
comparable and high [O in] A5007/H/3 ratios may still argue 
for separate ionizing sources, but it is possible that the nu- 
cleus associated with the blueshifted component is ionizing 
both velocity components. Chandra observations would clar- 
ify the situation considerably. To summarize, our imaging and 
spectroscopy data strongly support the binary AGN scenario 
for at least three of the four objects. 

3. DISCUSSION 

As listed in Table [T] the four new kpc- scale binary AGNs 
that we have pres ented here have optical colors typical of 
type 2 AGNs (e.g.. iKauffmann et al.ll2003h . All four objects 
show tidal features suggestive of mergers ( 32.11 Figure 
Among them, SDSSJ1 10851.04+065901.4 was detected (but 
unres olved) in FIRST w ith an integrated flux of 9.84 ± 0.13 
mJy dWhite et alJ Il997[) : it also has a candidate m atch in 
the IRAS Faint Source Catalog dMoshir et al.lll992h . IRAS 
Fl 1062+0715, with L IR ~ 1O 12 L , although the positional 
match is only marginal (~9.5 // offset from its SDSS position). 
SDSSJ1 13 126.08 - 020459.2, SDSSJ1 14642.47 + 5 1 1029.6, 
and SDSSJ133226.34+060627.4 were not detected in FIRST 
or in IRAS. At least three of the four galaxies appear to be 
neither star-bursting luminous IR galaxies (LIRGs) nor pas- 



sively evolving ellipticals. These properties are in contrast to 
the few existing samples of known kpc- scale binary AGNs, 
which tend to be found either in LIRGs which involve gas- 
rich mergers of disks (N GC 6240, Komossa et al. 2003; Mrk 
463, iBianchi et all 120081: Arp 299. IB alio et al.ll2004 . or in 
ellipticals a t the centers of galaxy clusters a nd power twin 
jets (3 C 75. IOwen et aDfl985l: PKS 2149-158. iGuidetti et al.1 
2008). Perhaps this reflects our differing selecting criteria, but 
we are hopeful that with our systematic approach we can be- 
gin to characterize the evolutionary state of kpc- scale binary 
AGNs in more detail. Further host galaxy studies on these 
four binary AGNs including their structural properties, stellar 
populations, and dynamics, as well as the remainder of our 
NIR imaging subsample will be presented elsewhere. 

Until recently, searches for kpc- scale binary AGNs have 
been more or less serendipitous. We have demonstrated with 
the observations presented in this paper that selecting spec- 
troscopic candidates based on double-peaked narrow lines 
and follow-up with high spatial resolution imaging and spec- 
troscopy is a promising technique to identify genuine binary 
AGNs in an efficient and systematic way. Out of the 43 ob- 
jects observed so far, we detected three, and perhaps four bi- 
nary AGNs (~ 7-10%). These four systems probe the rel- 
atively early stages of AGN pairing, when the two stellar 
bulges which host the two SMBHs are still separated enough 
to show two distinct components in IR/optical images. The 
~ 10% lower limit translates into a 0.1% kpc-scale binary 
AGN fraction out of the parent AGN sampl e (167 double - 
peaked objects selected from 14,756 AGNs iLiu et~ai1l2010l) . 
comparable to the binary qu asar fraction with projected sepa- 
rations of tens of kpc (|Hennawi et al.l2006t iMyers et ail 200 8 ; 
Henn awi et al.l2009|) . albeit with different systematic s and po- 
tential selection biases. We also have objects in our sam- 
ple in which the [O in] emission is spatially resolved but the 
NIR continuum is not, and in a few such cases the spatially- 
resolved spectrum shows clear signs of a velocity gradient 
across the slit, suggesting a rotation/outflow origin for the 
double-peaked profile. Thus it is premature to claim a binary 
AGN based on spatially-resolved [O in] emission alone. A 
detailed analysis on the fraction of such "false positives" and 
their properties will be presented in a future paper. 

We are carrying out NIR imaging of the rest of our double- 
peaked narrow-line sample using Magellan and Keck with 
AO. Together with follow-up spectroscopy, we expect to in- 
crease the number of known kpc-scale binary AGNs several 
fold. High-quality multi-band data will allow us to explore the 
properties of individual binary AGNs in great detail, and the 
increased statistics will provide important constraints on the 
merger hypothesis, AGN lifetimes, and the interplay between 



BINARY AGNs 



5 



SMBHs and their hosts. 
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